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Persistent, oscillatory charge and spin currents are shown to be driven by a two-component 
terahertz laser pulse in a one-dimensional mesoscopic ring with Rashba-Dresselhaus spin orbit inter- 
actions (SOI) linear in the electron momentum. The characteristic interference effects result from 
the opposite precession directions imposed on the electron spin by the two SOI couplings. The time 
dependence of the currents is obtained by solving numerically the equation of motion for the density 
operator, which is later employed in calculating statistical averages of quantum operators on few 
electron eigenstates. The parameterization of the problem is done in terms of the SOI coupling 
constants and of the phase difference between the two laser components. Our results indicate that 
the amplitude of the oscillations is controlled by the relative strength of the two SOI's, while their 
frequency is determined by the difference between the excitation energies of the electron states. 
Furthermore, the oscillations of the spin current acquire a beating pattern of higher frequency that 
we associate with the nutation of the electron spin between the quantization axes of the two SOI 
couplings. This phenomenon disappears at equal SOI strengths, whereby the opposite precessions 
occur with the same probability. 

PACS numbers: 73.23.Ra,71.70.Ej,72.25.Dc,73.21.Hb 



I. INTRODUCTION 

The reconsideration of the spin-orbit interaction (SOI) 
in the modern context of spintronic applications, started 
more than ten years ago, was based on the idea of the 
possible manipulation of the electron spin by means of 
an electric fieldji While definitive answers to this quest 
have yet to be reached, the fundamental interest in under- 
standing the effect of SOI on macroscopic phenomenol- 
ogy continues. Associated with the spatial confinement in 
two-dimensional quantum wells^ (Rashba-R) or with in- 
version asymmetry in the crystal (Dresselhaus-D)^, SOI 
is usually considered to be linear in the electron momen- 
tum which is coupled to the electron spin through a con- 
stant whose magnitude is amenable to outside control. 
The most important physical aspect associated with the 
R-D superposition is the effect on the electron spin, which 
is forced to precess in opposite directions along the quan- 
tization directions imposed by the two spin interactions. 
This phenomenology is embodied by the single particle 
Hamiltonian, 

a P 

HsO = -^{<^xPy - <yyPx) + -^{(^xPx ~ CFyPy) , (1) 

which introduces the Rashba and Dresselhaus interaction 
constants a and /3 which couple the electron momen- 
tum p = (pxiPy) with the electron spin represented by 



the Pauli spin matrices, (Jx,y The superposition of the 
Rashba and Dresselhaus terms generates particularly in- 
teresting situations when their strengths are equal, such 
as the cancellation of dephasing for the eigenstatc spinors 
and the ensuing ballistic spin transporlj^ or the formation 
of a persistent spin helix.— 

In this paper, we investigate the consequences of 
the Rashba-Dresselhaus superposition on spin and 
charge currents that are being induced in a quasi-one- 
dimcnsional ring by a terahertz laser pulse with a spatial 
asymmetry. This represents a well known method of cur- 
rent generation which exploits the left-right asymmetry 
of the electron states that arc excited on a time scale that 
is much shorter than the electron relaxation lifetime in a 
non-adiabatic process.— In general, by this method, as 
well as by applying magnetic fieldsf^i^ simultaneous spin 
and charge current generation occursJ^^— The indepen- 
dent realization of pure spin currents has been addressed 
in a number of theoretical proposals that considered hy- 
brid structuresi^ or a specific electron configuration (odd 
numbers of particles) More recently it was shown that 
a pure spin current can be created non-adiabatically in a 
ring with Rashba interaction using a radiation pulse with 
two dipolar components having a spatial dephasing angle 
0.— The physical mechanism for spin current generation 
relies on the interplay between the spin orbit coupling 
that rotates the electron spin around the ring and the 
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spatial asymmetry of the external excitation, which es- 
tablish conditions where the charge current disappears, 
while the spin current reaches a maximum or a minimum 
level. 

This work undertakes the analysis of spin and charge 
dynamics in the simultaneous presence of Rashba and 
Dresselhaus SOI, such that two opposite precession di- 
rections arc imposed on the electron spin. Our numeri- 
cal results arc obtained within an equation-of-motion al- 
gorithm for the particle-density operator, that is later 
involved in calculating the spin and charge currents as 
statistical averages of the corresponding quantum oper- 
ators on a few non-interacting electron eigenstates. The 
internal phase difference </> between the two laser com- 
ponents along with the coupling strengths of the R and 
D interactions are the main parameters of the problem. 
While the latter play a role in determining the magni- 
tude of the effects, the former is shown to influence the 
spatial distribution of the currents around the ring. 

The paper starts by discussing the spectrum of the 
equilibrium Hamiltonian, in Sec. whose eigenstates 
and eigenvalues are obtained within a direct diagonaliza- 
tion procedure for a small number of electrons. Then, in 
Sec, mil wc detail the equation-of-motion algorithm that 
allows the estimation of the density operator, followed by 
Sec. IIVI and fVl that present our findings. The results are 
summarized in Sec. IVII 



II. THE EQUILIBRIUM HAMILTONIAN 

Our analysis is based on a discrete model of the Hamil- 
tonian that relies on transforming the continuous, quasi- 
one dimensional ring of radius into a sequence of 
N sites (points) separated by an equivalent lattice con- 
stant a = 2'Kro/N. In polar coordinates a site of in- 
dex n = 1,...,A^ has an azimuthal angle 0„ = 2Trn/N, 
the angle difference between two consecutive sites being 
A9 = 2-k/N . A single particle state that corresponds to 
an electron of spin a located at point n, cr), is asso- 
ciated with the creation and annihilation operators c!^ ^ 
and Cn,a- The electrons encased in the ring are described 
by the Hamiltonian H^ng composed out of the kinetic 
energy Ho and the Rashba and Dresselhaus terms, Vr 
and Vd, 

Hrins = Ho + Vr + Vd , (2) 

where Vr and Vd are given by the first and the second 
term of Eq. ([T]), respectively. In the representation pro- 
vided by the single particle states described above, Hq is 
proportional to the hopping energy V = h^/2m*a^ (m* 
being the effective mass of the electron in the host semi- 
conductor material). 

Ho = 2F^c+^c„_^-V^^c+^c„+i^^-V^^c+^c„_i^^ , 

n,(T n,a n,cr 

(3) 



while the SOI terms Vr and Vd generate their own energy 
scales, Va = a/2a and Vp = /3/2a, 

VR = -iVc J2 ['^r{0n,n+l)]„^„,C+^C„+la' (4) 

and 

VD = -iVp [ct9(6I„,„+i)]*_^, c+^c„+i^, -hH.c. . (5) 

For simplicity, it is customary to introduce the azimuthal 
and radial spin matrices ag and written in terms of 
the angle 0„^n+i = {On + 6'ri+i)/2, 

<7r{0)= ax COS 6 + <7y sin , 

ag{9)= —axSinO + (TyCosO . (6) 

In the presence of only one type of SOI, say Rashba, the 
spectrum of the Hamiltonian, as well as the eigenvalues 
of the spin operator, can be obtained analytically. The 
single-particle eigenvectors are also eigenstates of L^, the 
z component of the angular momentum, and they are 

\^,^) = ^ye^^o,J cose \ 

with ; = 0,±1,- - • ,±(iV- l)/2,iV/2 (for N even). Cor- 
respondingly, the energy eigenvalues are given by 

i?.± = ^^^ + ^^^v/l + tan^20„, (8) 

where ei = 2V - 2V cos{2ttI / N) and tan26'„ = 
Va/ {V sv[\{tt /N)). The quantization direction of the spin 
operator is e2e„ = cos 20a.ez — sin 20Q,er , which is tilted at 
an angle 29a relatively to the ^-axis due to the presence 
of the of the SOL— The eigenstates are spin degenerated 
and the energies Ei± correspond to spin eigenvalues of 
±h/2 along e2ea- In turn, the Dresselhaus interaction 
alone determines a similar energy spectrum and eigen- 
states, but defines another preferential spin orientation 
direction, 620^ ~ cos26'/je2 -|-sin26'^eg where = — e_e 
(e_e is azimuthal direction for the point of angle —6 of 
the ring). 

When both R and D interactions are present, the di- 
agonalization of the Hamiltonian ([2]) is possible only by 
numerical methods and the results are strongly affected 
by the relative values of the two coupling constants a and 
p. The energy spectrum of the equilibrium Hamiltonian, 
calculated for an even number of sites, N = 20, is shown 
in Fig. [T] as a function of Va , for different values of Vp . 
All states remain spin degenerate, but other degeneracies 
are lifted. When Vg = 0, as in Fig.[T][a), the spectrum 
is four fold degenerate for those values of Va that allow 
the equality Ei^ — Eii_ for two different quantum num- 
bers I and I' . Such degeneracies are lifted when Vp ^ Q 
and comparable to Va, as shown in Fig.[IJb), where we 
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III. THE NUMERICAL ALGORITHM 

The perturbation Hamiltonian i/puisci which repre- 
sents a laser pulse with two dipolar components of fre- 
quencies uji and UJ2, out of phase with each other by an 
angle 0, is acting on the few rie = 6 electrons that occupy 
the lowest energy eigenstates. In the discrete represen- 
tation of the ring, we write J?puiso(i) = 
where the time dependence Hg{t) at the point n of the 
ring is: 

Hg {t) = Ae~^*[sm{uJit) cos 0„ + sin(a;2t) cos(0„ + cj))]. 

(11) 

The pulse of amplitude A is applied at time t ~ and 
lasts for a time i/ ~ F^^. The density operator p{t) 
satisfies the differential Liouvillc equation, 



ihp{t) = [H + Hp^u,{t), p{t)] , 



(12) 



whose solutions are obtained numerically, subject to the 
initial conditioi*^ 



(13) 



1=1 



FIG. 1. (Color online) The energy spectrum of a ring with 
SOI versus Rashba energy Va- The Dresselhaus energy is 
Va = in (a) and Va — 0.05 in (b). All energies are given 
in units of V = /2m*a^ . The spectrum has 40 eigenstates, 
corresponding to 20 sites x 2 spin directions. All eigenstates 
are spin-dengenerated. 



notice that small energy gaps appear between previously 
intersecting levels. These gaps give rise to low frequency 
oscillations of spin and charge currents. 

The azimuthal velocity operator vg — r^dO/dt = 
(iro I h) [H, 9] is evaluated from Eq. ^ , 



n.(T 



+ 2^ E ^r{OnY^"' cl„C„ 



+2^ E ^iOuT^"' cl„c.^.„, 



(9) 



leading to the charge and spin current operators: 

/ = -evg , 



(10) 



In Eq. (|10p . a I, is a spin matrix corresponding to the 
direction e^^ which can be any of the Cartesian vectors 
^x, Bj,, or e2e„, 62^^. 



where only the linear superposition of the lowest rig oc- 
cupied states is considered. For any t > Eq. 
generates pit) through the Crank-Nicholson finite differ- 
ence method^ii^ with small time steps 5t <^ F~^. The 
expected value of any observable O of the many-body 



p{t)0 



where O is the 



system is then calculated as Tr 

single-particle quantum operator associated to that ob- 
servable. After the perturbation ceases to act, at time 
t > tf, the system remains in an excited state of constant 
energy. The time evolution of the expectation values of 
the system observables is determined by their commuta- 
tion relation with the Hamiltonian. 

In the following calculations the radius of the ring is 
tq = 14 nm, which for = 20 sites generates a lattice 
constant a w 4.4 nm. The effective electron mass in the 
ring is material dependent, m* = O.OBTme for a GaAs or 
m* = 0.023TOe for InAs. (?7ie is the free electron mass.) 
Table 1 summarizes the various parameters involved in 
the calculation that pertain to the applied laser pulse 
(time, frequencies, pulse amplitude A) and those that 
pertain to the physical observables of the system (energy, 
spin, velocity, charge and spin current Ic and Is). We also 
show the correspondence between the SOI parameters Vr 
and Vd and Rashba and Dresselhaus constants a and f3 
in the two materials. 

In the numerical calculation we use the pulse frequen- 
cies wi = O.G963y/;i and uj2 = 0.276V/h, while the atten- 
uation factor and the amplitude are chosen as F = 4a;i 
and A = 2.3V, respectively. The two frequencies are 
close to the first two Bohr frequencies of the ID ring 
calculated in the absence of SOI, u}2i ~ 0.0978y/?i and 
UJ32 ~ Q.284V/h. A Bohr frequency utij is given by the 
energy difference (i?2i — E2j)/fi, with i?2i and E2j being 
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Physical units 



Parameter 



Unit 



GaAs 



InAs 



Energy 
Time 
Frequency 
Velocity 
Charge current 
Spin current 



V = ft7(2m*a^) 
hjV 
V/h 
aV/h 
eaV/h 
aV/2 



29.4 meV 
0.022 ps 
44.6 THz 



85.6 meV 
0.0076 ps 
130.0 THz 



196-10^2 nm/s 572-10^2 nm/s 
31.5 /iAnm 91.6 /^Anm 
64.6 meVnm 188 meVnm 



TABLE I. The specific units for energy, time, frequency, 
charge and spin currents, for GaAs (m* = 0.067me) and for 
InAs (m* = 0.023me) quantum rings, used in the numerical 
calculations, e is the electron charge. The discretization con- 
stant is a=4.4 nm and is calculated for a ring of radius ro = 14 
nm with N=20 discrete points. The charge current is defined 
as Ic = ev and spin current is Is = sv, v being the velocity 
and s the average spin angular momentum in units of h. 



the 2i*'' and 2f^ eigenvalues of the quantum ring. The 
factor 2 takes into account the spin degeneracy. The ex- 
ternal pulse last about = lO/F and is equal to 2'ih/V 
dimensionless time units. 

According to Table U for a GaAs quantum ring, the 
above values correspond to energies hwi = 2.83 meV, 
huj2 = 8.11 meV, the attenuation factor F = 4a;i, the 
amplitude A = 67.7 meV and the external pulse duration 
of about = 0.2 ps. If an InAs quantum ring is consid- 
ered, the numerical values associated with the laser pulse 
arc hoji = 8.23 mcV, fiaj2 = 23.6 meV, A = 197 meV and 
tf = 0.068 ps. 

The strength of the SOI is chosen by using the pa- 
rameters Va and Va in interval [0, 0.08^], where V is the 
energy unit shown in Table HI A value Va = 0.05T^ corre- 
sponds to a Rashba coupling constant a ~ 12.9 meVnm 
in a GaAs quantum ring or a = 37.6 mcVnm in InAs as 
shown in Table [Til 

Values of Rashba and Dresselhaus parameters 



Definition 



GaAs 



InAs 



a = 2aVc a=258.4 meVnm-Vc/V q=752.7 meVnm-Vc/V 
l3 = 2aVp /?=258.4 meVnm-Va/V =752.7 me Vnm-Va/V 

TABLE II. The Rashba and Dresselhaus parameters a and 
/3 for GaAs and InAs quantum rings are expressed in terms of 
the Rashba and Dresselhaus energies Va and Vg, respectively, 
in units of V . The calculations are done for a ring of radius 
ro = 14 nm with N=20 discrete points. 



IV. CHARGE AND SPIN CURRENTS DRIVEN 
BY A SYMMETRIC PULSE 

We begin our investigation by analyzing the behavior 
of charge and spin currents in the presence of an R-D SOI 
when the ring is excited by a spatially symmetric laser 



pulse, i. e. (/> = in Eq. (|TT|) . We note that at time t ~ 0, 
the charge current is zero for all strengths of SOI, Va or 
Vp. Under the effect of the perturbation, left-right elec- 
tron state imbalance occurs, leading to the establishment 
of an oscillatory current, as shown in Fig. [2l with period 
T of thousands of time units h/V. There, the inten- 
sity Ic{t) is plotted for various Rashba couplings when 
the Dresselhaus coupling assumes two different values, 
Vp = OmV in Fig. [ila) and Fig. ^b), and Vp = 0.051/ 
in Fig. ^c) and Fig. [IJd). The period and amplitude 
of the current depend on the SOI parameters, reaching 
a maximum when Va = Vp. The small oscillations of 
amplitude A'l^ < 0.012 eaV/h and period T' « 16 h/V 
noticeable in Figs.[2{c) and[2][d) are further investigated 
in the discussion of the spin currents. 

The low frequency of the charge current O = 1/T is 
given by the Bohr frequency corresponding to the occu- 
pied states with highest energies in the unperturbed ring, 
Eq and E4, such that fl — 0132/277. We verified this re- 
sult for all the SOI strengths used in these calculations. 
For example, for Va = 0.05 and Vg = 0.03 the first few 
(relevant) energy levels of the quantum ring are, in the 
increasing order, Ei ^ E2 = -0.003383, E3 = E^ = 
0.09131, E^ = Eq^ 0.09805, Ej = Es = 0.3735, Eg = 
EiQ = 0.3850, Ell = E12 = 0.8146. We thus obtain from 
the energy spectrum T = 27r/a;32 = 27r/ {Eq — E4) = 932 
time units, whereas directly from the numerical results 
of the time dependent current we get T = 948. This 
means that after the perturbation had ceased, the ring 
remains in an excited state which is a superposition of 
states with partial population of the energy levels Eg and 
E4. The frequency of the oscillations, vs. the Rashba 
SOI Va e [0, 0.08V], is plotted in Fig.O for two values of 
the Dresselhaus SOI. Both the results obtained directly 
from the time dependent charge current, and from the 
Bohr frequency UJ32, are shown. The agreement is al- 
most perfect. We find a non-monotonic behavior that 
is highly dependent on the relative strengths of the two 
SOI's, such that for a fixed Vg, the minimum frequency 
occurs for equal strengths Va = Vp. In this case addi- 
tional degeneracies exist in the energy spectrum^ and 
therefore when Va Vp all energy gaps tend to shrink. 
(With our parameters, for equal SOI strengths, we obtain 
Ei5 = EiQ = Eij = Eis.) 

The fast oscillations of the charge current, which can 
be seen well in Figs. [DJc-d) can also be explained by the 
energy spectrum. They are created by transitions be- 
tween states separated by relatively large energies, like 
Eq, E4, E2 and Eg, or Eiq. In fact the period of the fast 
oscillations corresponds to the largest energy gap in this 
series, T' = 2n/{Eio - E2) ~ 16 h/V. 

Under the same circumstances we study the time vari- 
ation of the spin current associated to the spin projec- 
tion along the z-direction, (t) , for various values of the 
SOI interactions. In this instance, the oscillatory behav- 
ior associated with charge displacement is complicated 
by a a beating pattern with nodal points of zero oscil- 
lation, as seen in Figs. H^a) and (c). Zooming in near 
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FIG. 2. (Color online) Charge currents Ic{t) in the presence 
of both Rashba and Dresselhaus SOI. V/j = 0.03 in (a) and 
(b) and Vp — 0.05 in (c) and (d). Va is written in the figures. 
The minimum frequency Qmin ~ 1/Tmaa; and the maximum 
amplitude of the oscillations in time, Ale, are obtained for 
Va^Vfi. Qmin = 0.00028 for Vc = Va = 0.03V and n^i„ = 
0.0008 for Va — V/i — 0.05V. The ring and pulse parameters 
are described in the text. Ic is expressed in eaV/h and time 
in h/V units. 
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FIG. 3. (Color online) The frequency of the long oscillations 
of the charge current are plotted as functions of the Rashba 
energy Va for two values of the Dresselhaus coupling l/jj = 
0.03V and V/j — 0.05V respectively. The continuous lines are 
obtained from the time dependent output data shown in Fig. 
(2] The dotted lines show the results for the Bohr frequency 
i^32/27r = {Eq — E4) / {2n) (see text). The minimum frequency 
is obtained for Va = Vp. The frequency unit is V/h and Va is 
in units of V. 



the nodal points, wc obtain Fig. |4l^b) and|4|d), where it 
is readily observable that, as it passes through a nodal 
point, the spin current phase changes by tt. As before, 
the amplitude of the oscillations depends on the two SOI 
interactions and their ratio, while the period of the fast 
oscillations T' seems quite insensitive to this aspect. The 
period of the fast oscillations is T' « 8 h/V or slightly 
more for all examples shown in Figs. HJa-b), while the 
amplitude varies from zero, at nodal points, to about 
0.1 aV/2. The period T' of the spin current is actually 
half of the period of the fast oscillations of the charge 
current (also denoted by T' in the analysis of that cur- 
rent, and found there « 16 time units). The reason is 
that the spin precesses relatively to two principal axes, 
independently. The axes are tilted at different angles 26*0, 
and 29 due to the Rashba and Dresselhaus SOI, respec- 
tively. The frequency of the spin current is thus roughly 
the double of that for the charge current, because of the 
two spin modes For Va = Vg the spin oscillations can- 
cel each other and the spin current vanishes, as seen in 
Figs.mja) and (c). For V^^Vp the spin current actually 
has an aperiodic time evolution, observable in the same 
mentioned figures. So we cannot really speak about a rig- 
orous periodicity of the spin current, and consequently of 
the charge currents as well. 

We interpret the beating pattern of the spin current 
oscillations by a nutation motion of the electron spin 
between the two quantization directions ^26^ ^-iid e2e^ 
imposed by the SOI couplings that force the electron 
spin precession in opposite directions. When only one 
type of the SOI, say Rashba, is present, e2e„ becomes 
a good quantization axis for spin. The corresponding 
spin current 128^ commutes with the equilibrium Hamil- 
tonian and consequently, remains constant in time as 
a conserved observable. When the SOI parameters are 
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FIG. 4. (Color online) The time variation of the spin current 
Iz{t) in the presence of Rashba and Dresselhaus SOI. The 
Dresselhaus energy is Va = 0.03 in (a) and (b). The Rashba 
energy is Vc = 0.01, 0.02, 0.03 in (a) and K = 0.03, 0.04, 0.05 
in (b). For Va = Vis the spin current is zero. For Va 7^ Vg 
there are spin current oscillations. Figures (c) and (d) present 
a magnified view of (a) and (b) respectively. (Spin current 
unit is aV/2 and time unit is h/V .) 



not equal, the principal spin axes have different angles 
26*0, 7^ 2013 and the electron spin executes rapid oscilla- 
tions between them. 

Another perspective on this phenomenon is revealed 
by plots of the time oscillation of the SOI energy, Eji 
and Ed, defined as the expected values of the potential 
energies Vr and Yd given by Eqs. ^MSi- Figs. [Sfa) and 
mb) display ER{t) and ED{t) for different values of SOI, 
while ininija) andlH^b) the same value Va = Vg is used. It 
appears that an exchange between Rashba and Dressel- 
haus potential energies occurs with the same oscillation 
period as the Iz beat oscillations. The orbital motion 
of the electron is thus accompanied by vibrations of the 
spin orientation between these two directions that give 
rise to out of phase oscillations of Rashba and Dressel- 
haus energy. At equal SOI parameters, Va = Vg, the two 
potential energies reach the same value and oscillate in- 
phase. Since the precession of the electron spin around 
the two preferential directions or the Rashba and Dres- 
selhaus SOI's occurs with equal amplitudes in opposite 
directions, the spin current vanishes. 



V. CHARGE AND SPIN CURRENTS DRIVEN 
BY AN ASYMMETRIC PULSE 

The time evolution of the electronic states in the quan- 
tum ring excited by an asymmetric pulse is asymmetric 
is discussed next. This involves the full form of Eq. PT|) . 
with 7^ 0. In this case, we are interested in establishing 
the geometric effect that the dephasing angle 0, varied 
in fixed steps = 0, 7r/6, 27r/6, • • • , ll7r/6, has on the 
charge and spin currents. To better focus on this aspect, 
the Rashba and Dresselhaus interactions are fixed at val- 
ues Va = 0.02y and Vg = 0.03T^, respectively. In our 
analysis, two dynamic regimes are considered, delimited 
by the lifetime of the perturbation, tf, one occurring for 
t < tf, and the other one for t > tf when the persistent 
oscillatory behavior is established. 

In Fig. [7{a) the charge current Ic{t,4>) is shown to 
evolve under the effect of the two-component laser pulse 
between < t < tf, from zero at t = to non-zero val- 
ues at t = i/, for several angles (j>. The largest currents 
occur for dephasing angles (j) = An/G and 8tt/6, whereas 
the minima occur for = and tt. A similar pattern is 
seen in Fig.[7Jb), where the dependence of the instanta- 
neous charge current at a given time /c(t,0) is depicted 
for various angles 0. A maximum variation is noticeable 
at = 7r/2 and 3tt/2 only for t = 5h/V and t = 8h/V, 
while for larger times, such as at f = 20h/V for instance, 
the maximum variation of the induced current |/c(OI 
realized for dephasing angles = Air/ 6 and Stt/G as al- 
ready seen in Fig. [71[a). 

After the radiation pulse vanishes, i. e. for t > 
t/ ~ 25h/V, the charge current oscillates with a period 
T ~ 2783h/V, as shown in Fig. [TJc), where now the 
time covers one complete period. The amplitude of the 
oscillations is strongly affected by the value of the de- 
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FIG. 5. (Color online) The out of phase oscillations of the 
Rashba and Dresselhaus potential energies, ER{t) and Eoit). 
Vjs = O.OSy and Va = 0.021/, (a) for the complete time inter- 
val, and (b) zoomed into a small time interval. The energy 
unit is V and time unit is h/V . 
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phasing angle. For the dephasing angles = 0, 7r/6 and 
ll7r/6 the time dependence of the charge current has the 
shape of the function — sin(27rt/T), and will be called 
here type A oscillations. For the angles (f) = 37r/6, 47r/6, 
87r/6, and 97r/6 the oscillations look like sin(27rt/T) and 
will be called type B oscillations. Type B are shifted 
from type A by half a period, T/2. Type A oscilla- 
tions of much smaller amplitudes are also seen in Fig. 
[Tjc) for = 67r/6, and also for (j) = 27r/6 and IOtt/G. 
The latter two are actually asymmetric and in transi- 
tion from type A to type B. For </> = 57r/6 and 77r/6 
again small amplitudes are obtained, this time for oscil- 
lations of type B. The amplitude of the current, defined 
as A/c = |/c(T/4) - /c(3T/4)|, changes with 4> as shown 
in Fig. [71[d). As the angle is varied, the two types of 
oscillations alternate. This result defines four critical an- 
gles (f) — 0ci,c2.c3,c4 that correspond to a crossover from 
type A to type B. Their numerical values obtained by the 
interpolation of the presented curves are: O.SStt, O.SStt, 
I.IStt, and 1.657r, respectively. At the A to B crossover, 
1. e. at the critical angles 0ci,...,c4, the oscillations of 
the charge current vanish, as shown in Fig. [71[d). This 
means that, by tuning the value of the dephasing an- 
gle of the external pulse such that 4> = 4'ci,...,ci the ring 
can be excited into a final state that supports a non-zero 
charge current with zero amplitude, although the cur- 
rent operator Ic is a nonconservative observable in the 
presence of both Rashba and Dresselhaus interactions. 
The time average of the charge oscillations, defined as 
{Qt = (/c(T/4)-f /e(3r/4))/2, is also shown in Fig.EId) 
as function of 4>. We note that, for type A oscillations, 
{Ic)t = for angle = 0m = and tt, and acquires 
its maximum and minimum values of ±0.4 at the angle 
= 0J\/ = O.Gtt and = (I>-m = 1.47r, for B type oscilla- 
tions. This suggests that, by tuning the dephasing angle 
4>, both the type of charge current oscillations, as well as 
their time average can be modified. 

In a parallel analysis, we present the time evolution of 
the spin current in Fig.[51 At time t = 0, the ring is found 
in the ground state, which has a permanent spin current 
4(i = 0) = 0.01760T/a/2. Along the SOI principal axes 
the spin currents are l2e^{t = 0) = 0.017751/a/2 and 
l20p{t = 0) = 0.01793l/a/2, respectively. After the onset 
of the laser pulse, for t < tf, the spin current changes 
non-adiabatically, as displayed in Fig. [DJa) for different 
dephasing angles (j). In this time interval the magnitude 
of Iz{t) decreases in time, the variation 1^(0) — Iz{t) being 
a function of </>. Snapshots of Iz{t,(t)) recorded at times 
t = 0, ^h/V, %h/V are plotted against the dephasing 
angle (f) in Fig.[51[b). We note that the smallest variation 
-^2(0) — Iz{t) occurs for = and the largest for = 27r. 

After the radiation pulse vanishes, i. e. for t > tf, 
the quantum ring enters the oscillatory regime. Similar 
to the charge current, the spin current has an oscilla- 
tory behavior with period T ~ 2783h/V shown in Fig. 
[51[c) where we plot Iz (t) . The amplitude AIz varies from 
—0.6Va/2 to 0.6Va/2 when the angle (j) of the external 
pulse is changed from 2tt/6 to IOtt/G, as seen in Fig. 



FIG. 6. (Color online) The in phase oscillations of the Rashba 
and Dresselhaus potential energies, Enit) and Eoit), for 
Vi3 — Va — 0.03V, (a) for the complete time interval, and 
(b) zoomed into a small time interval. The energy unit is V 
and time unit is h/V. 
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FIG. 7. (Color online) (a) The charge current Ic{t, (f>) as func- 
tion of time during the radiation pulse for t £ [0, 33], for dif- 
ferent angles = j'tt/B between the two dipoles. The integer 
j = 1, 2, 11 is indicated near each curve. In (b) we show the 
charge current as a function of at several moments shown 
in the graph, t = 0, 5, 8 and 20. (c) The time evolution of the 
current over a full period of the oscillation T « 2783. The 
angles are the same as in (a), (d) The average charge cur- 
rent (/c)t and the amplitude Ale are shown versus the angle 
0. The time average is {Ic)t = [Ic{T/A) + Ic{3T/i)]/2 and the 
amplitude is Ale = /c(r/4) - 7c(3r/4). The intervals of type 
A and type B oscillations are indicated. The SOI parameters 
are Va = 0.02V and V/s — 0.03V. Time physical units are 
h/V for the time and eaV/h for the current. 
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FIG. 8. (Color online) (a) The spin current Iz{t,(j)) during 
the pulse, for t £ [0, 33] in (a) for the same series of angles 
between the dipoles as in Fig. [Tja). In (b) we plot the spin 
current vs. the angle at several moments, t=0, 5 and 8. 
(c) The spin current along the z direction within a complete 
period T. The spin currents have long oscillations with period 
T ~ 2783 and short oscillations with T' ~ 8. The integers 
j — 1,2, ...,11 near the curve correspond to the angles — 
jn/6. In (d) we summarize the results by showing the time 
average of the spin current {I.)t = {I^{T/2) + J^(r))/2 and 
its amplitude AI^ = Iz{T/2) — Iz{T) versus the dephasing 
angle 4> ^-nd we indicate oscillations of type A and type B. 
The SOI parameters are Va = 0.02V and Va = 0.03V. Time 
unit is h/V and the unit for the spin current is Va/2. 



[5]Jc). In addition to the big oscillation w\X\v period T 
and large amplitude AI^ € [—0.6 : 0.6] Va/2, the spin 
current presents an overlapping pattern of small oscilla- 
tions, with smaller period T' ~ 8h/V and smaller ampli- 
tude A'lz — 0.03Va/2. A similar analysis can be done 
for the spin currents in the directions 28^ and 29/3 (not 



shown in the figures). The small and fast beating oscil- 
lation have a lower amplitude than for Iz . In the present 
example the oscillations of 129^ are slightly larger than 
those of l20„ because Vg > Vq . 

The analysis of Fig. [SJc) indicates that there are two 
type of spin current oscillations. For the angle cf) = 
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(a) Charge current l(,(t,(j)) 




(b) Spin current ls(t,(^) 




FIG. 9. (Color online) By varying the dephasing angle (j) of 
the external pulse different types of charge and spin current 
oscillations are induced in the quantum ring, (a) The charge 
current has oscillations in time of type A or B. For type A 
the average current vanishes for the angles (f) between the 
dipolar components of the radiation pulse marked with m. 
Type B oscillations occur for (j) £ [4>ci, <j>c2] and (j> G [(jjcs, (/>c4]. 
In this case the average current has maximum positive and 
minimum negative amplitudes for the angles marked with M 
and -M. At the critical points <^ci- -4 the charge current is 
constant in time, (b) The oscillations of the spin current 
are also classified as type A and B. For type A the averages 
current has minimum and negative values at the point —Mi 
and — M2. For type B it has maximum positive values at the 
points M2 and Mi. The two types of oscillations, A and B, 
for both Jc and I3, are in antiphase, respectively. 



27r/6, 37r/6, 7r/6 and 7tt/6, the analytic behavior for the 
time dependence of spin current is fitted by a function 
cos(27rt/T) that defines type A oscillations. For values of 
the dephasing angle (/> = 57r/6, IItt/G, Qtt/G and IOtt/G , 
the analytic representation is — cos(27rt/T). This defines 
type B oscillations that arc T/2 dephased from type A. 
For certain values of (j> the amplitude of the broad oscilla- 
tions of Iz , ^Iz 7 become comparable with the amplitude 
of the fast and small oscillations of amplitude A'/^ . This 
happens for = 0, tt, 27r/3 and 47r/3 in Fig. IH^c). In 
the later case the large oscillations with long period T 
vanish and the remaining feature is the beating pattern 
with nodal points already shown in Fig. 21 



The amplitude of the big oscillation, A/^ = Iz{T/2) — 
Iz{T), is found to be negative for type A and positive 
for type B oscillations, which is shown in Fig. |Hl[d). By 
varying (j) the spin current switches between A and B type 
oscillations at four values of the dephasing angle, denoted 
as (t>si,...,si- As shown in Fig.[51Jd) the critical values are 0, 
27r/3, TT, and 47r/3 (j ~ 0, 4, 6, 8). The time average of the 
spin current, {Iz)t = {Iz{T/2)+Iz{T))/2, is also plotted 
in Fig.[8l^c). When </> changes {Iz)t varies between locally 
minimum and negative value for </> = (p-Mi = 27r/6 and 
for (f) = 4>-Ar2 = Ttt/G reached within type A oscillations, 
and locally maximum positive values for </) = (j}M2 = 57r/6 
and (j) = 4>Mi = 107r/6 within type B oscillations. These 
results indicate that the type of oscillations performed by 
the spin current induced by the two-component radiation 
pulse, as well as the average spin currents, can be selected 
by changing the dephasing angle 0. 

The behavior of the spin current along the directions 
e26(„ and 626/^ is qualitatively similar to the spin current 
Iz {t) , having the same periodicity and beating structure. 
The amplitude of the fast oscillations is however smaller 
than for Iz{t). 

In deriving these results we ignored the inhomogene- 
ity of the electron distribution around the ring, which 
is known to appear in the simultaneous presence of 
the Rashba and Dresselhaus terms J^i^^ Underlying this 
choice is the fact that the charge deformation in a realis- 
tic ring is small, and even questionable for a narrow two- 
dimensional ring, unless placed in an external magnetic 
field»i^ Moreover, when the electron-electron interaction 
is considered for more than two electrons, the charge fluc- 
tuation is flattened out due to screening For our one- 
dimensional ring model the electron density has minima 
at polar angles 7r/4 and 57r/4 and maxima at 37r/4 and 
77r/4. Therefore, since the circular symmetry is intrinsi- 
cally broken in the ground state, a radiation pulse with 
only one dipolar component would, in principle, be suf- 
ficient to induce persistent oscillations of the charge and 
spin currents 



VI. SUMMARY AND CONCLUSIONS 

We investigate the interference effect generated by the 
simultaneous presence of the Rashba and Dresselhaus 
spin-orbit interactions on charge and spin currents in- 
duced non-adiabatically in a quasi-one-dimensional ring 
by a two-component radiation (laser) pulse. Our nu- 
merical results are obtained for a system of few non- 
interacting electrons through a direct calculation that 
involves the exact, time-dependent solution of the den- 
sity operator. The main finding of this work is that the 
oscillatory behavior of the charge and spin currents is re- 
alized at a frequency equal to the difference between two 
excited energy states (Bohr frequencies). 

By varying the dephasing angle between the two 
dipoles of the external pulse, different types of charge and 
spin current oscillations are induced in the ring. The gen- 
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eral features are summarized in Fig. Ufa) and Fig. [9jb), 
respectively. By changing the internal dephasing angle 
4> in the two intervals [4>ci, 4>ci\ and {4)c2, </>c3] the oscilla- 
tions of the charge current in time are qualitatively like 
— sin(27rt/T), i. e. the oscillations of type A indicated in 
Fig. |9l[a). The average current is zero (i. e minimum) at 
the angles indicated by the letter "m" in Fig.|9l[a), which 
are = and ip = tt. For S [4>ci, 4>c2] and 4> S [0c3, 4'ci\ 
the charge current oscillates in time (qualitatively) like 
the function sin(27rt/T), which are oscillations of type 
B, phase shifted with T/2 relatively to the ones of type 
A. The time average has positive maximum and negative 
minimum values at the angles marked with M and —M 
in Fig. [9l[a). When the angle </> has the critical values 
4'ci,...,ci a crossover between the two types of oscillations 
occurs, and the charge current is constant in time. 

The spin current has also such oscillations in time. 
When the angle (j) is in the two intervals [0si,0s2] and 
[0s3 , (psAl the induced spin current oscillates in time as 
cos(27ri/T), indicated as oscillations of type A in Fig. 
ini^b). The time average of the spin current has the min- 
imum (negative) values at the points marked as —Mi 
and — M2. For G [0s2,</'s3] and G [(/)s4,0si] the spin 
current oscillate in time as — cos(27rt/T), marked as os- 
cillations of type B in Fig.[9lb), and dephascd with T/2 
relatively to type A. Their maximum positive values oc- 
cur at the angles marked with M2 and Mi. These are 
wide oscillations with long period T and large amplitude. 



In addition, the spin current has also tiny oscillations 
with a much smaller period T' and amplitude A'. These 
high-frequency oscillations are caused by the nutation of 
the electron spin between the spin axes imposed by the R 
and D couplings. When the angle (p is close to the criti- 
cal points 4'si,...,sA a crossover occurs between type A and 
type B oscillations (or vice versa) of the spin current. For 
these values of the angle (p the amplitude of the wide os- 
cillations of Iz (A) decreases and becomes comparable to 
the amplitude of fast oscillations (A'). The spin current 
oscillations reduce to a beating pattern, while the wide 
oscillation with long period T vanishes. 

After the original excitation disappears the system sus- 
tains two persistent types of oscillations of the charge and 
spin currents, which we classified as type A and type B. 
The two types of oscillations arc in antiphase. Transi- 
tions between these modes can be controlled by varying 
the dephasing angle between the two components of the 
radiation pulse, an idea with potential applications in the 
spin-based information technology. 
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